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ABSTRACT 

We discuss recent observations of accreting binary pulsars with the all- 
sky BATSE instrument on the Compton Gamma Ray Observatory. BATSE 
has detected and studied nearly half of the known accreting pulsar systems. 
Continuous timing studies over a two-year period have yielded accurate orbital 
parameters for 9 of these systems, as well as new insights into long-term accretion 
torque histories. 


1. INTRODUCTION 


There are over 30 known accreting pulsar systems. These contain a rotat- 
ing high- magnetic-field (B lO^G) neutron star in orbit with a stellar com- 
panion which transfers mass to the neutron star via Roche-lobe overflow or a 
stellar wind. The gravitational energy released by mass accretion yields thermal 
and non-thermal radiation, most prominently at X-ray and gamma-ray ener- 
gies. Neutron stars in accreting pulsar systems are often called. “X-ray pulsars” 
or “accretion- powered” pulsars. An earlier comprehensive review of accreting 
pulsars is that of Nagase (1989). 

In this paper, we discuss recent observational results on 14 of the accret- 
ing pulsar systems, using data from the Burst and Transient Source Experiment 
(BATSE) on the Compton Gamma Ray Observatory ( GRO ). We list these sys- 
tems in Table 1, together with current information on position, pulse period 
( P m ), orbital period (P or &), and stellar type of the companion. Figure 1 
indicates the location of the BATSE-studied pulsars and other pulsars in the 
( Pspin — Porb) plane, commonly known as the "C'orbet-diagram” (see Corbet 
1986; Stella, White, and Rosner 1986; and Waters and van Kerkwijk 1989). The 
three types of high-mass ( > 3 M@) systems (wind-fed, disk-fed, and Be) populate 
relatively distinct regions of the diagram. The continuous BATSE spin-period 
measurements allow tests of models of the long-term spin evolution of pulsars 
in all three classes of systems. 
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TABLE 1 

Accreting Pulsars Observed with BATSE as of 1993 December 
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For epoch TJD 8500. 

Suggested companion (Glass & Feast 1973; Davidsen et al. 1977; Chakrabarty et al. 1994) 

Companion not yet identified, but spectral type inferred from orbit and eclipse measurements (Chakrabarty et al. 1993) 
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Figure 1. Corbet Diagram. Wind-fed systems are indicated by squares, high- 
mass disk-fed systems by triangles, Be-binaries by diamonds, low-mass systems 
by crosses, and systems with an unidentified stellar companion by asterisks 


2. OBSERVATIONS 

Observations of accreting binary pulsars have been carried out by BATSE 
since the launch of GRO in April 1991. The observations used the eight large 
area detectors (LADs) of the BATSE instrument, which have overlapping fields 
of view covering a total of 47 r steradian. Each detector has an effective area 
of about 1500 cm 2 at 40 keV and an energy resolution of about 35% FWHM 
(Fishman et al. 1989). The LADs are sensitive from about 20 keV to 2 MeV 
and pulsar timing observations are primarily carried out from 20 keV to 60 keV . 
Several data types were used for pulsar analysis: 

DISCLA data: Count rate samples of all 8 detectors at 1.024 s intervals 
in 4 energy channels. 

CONT data: Count rate samples of all 8 detectors at 2.048 s intervals in 
16 energy channels. 

PSR data: Count rate samples folded into 64 phase bins with a pro- 
grammable folding period and 16 energy channels (4 for periods shorter 
than about 20 ms). 

DISCLA and CONT data have been used for timing and spectral studies of 
pulsars with pulse periods greater than about 2 s and 4 s respectively. PSR 




data have been used for study of fast accreting pulsars, in particular Her X-l. 
The DISCLA and CONT data provide continuous flux information on all pulsars 
not occulted by the earth, while the PSR data have typically been programmed 
to study specific pulsars with known pulse periods. 

To detect a pulsar, the count rates from detectors with significant pro- 
jected area towards the pulsar are weighted and summed. Systematic orbital 
background is subtracted either by filtering or by use of a background model. 
The power spectrum of the resulting residual counting rate is calculated and 
searched for significant peaks. Once a pulsar is detected, timing measurements 
are carried out in the usual fashion by constructing a pulse-phase or pulse-arrival 
time model which accounts for barycenter corrections, the orbit of the pulsar, 
and torque-induced spin-frequency changes. The best-fit model is used to de- 
termine the orbital parameters of the pulsar system. 

The measured sensitivity of the BATSE detectors for pulsed flux detection 
is shown in Figure 2. Searches up to this time (Jan 1994) have used 1-4 day 
observation periods. The sensitivity is affected by earth-orbital noise for periods 
greater than about 200 s. The sensitivities for individual CONT energy channels 
are indicated in Figure 2, as well as sensitivities obtained by combining several 
energy channels. Typically, the lowest three CONT channels shown in Figure 
2 are used for pulsar detection. Also shown for comparison is the approximate 
pulsed Crab flux; the total Crab flux is a factor of about 5-10 higher depending 
on energy. 


10.00 I 


o 

CD 

E 

u 

> 

CD 


Crab RMS Pulsed Flux 

1.00 b ; 


X 

=3 


CD 

in 

13 

CL 

c n 

DC 

X 


o,iob— e- 


— e — e- 


-o — 


3a Sensitivity 

100000 Second Observation 


0.01 


100 


1000 


Energy (keV) 


Figure 2. BATSE Pulsed-Flux Sensitivity. 




3. RESULTS 


Dynamical Studies 

BATSE has determined the orbits of 0A0 1657-415 and GS 0834-430 and 
has provided improved or additional orbital solutions for < other systems: 3 
Be-binary systems (EXO 2030+375 A0535+26, and 4U 0115+63), 3 high-mass 
supergiant systems (Cen X-3, Vela X-l, and 4U 1538-52), and one low-mass 
system (Her X-l). In addition, one new accreting pulsar has been discovered, 
GRO J 1008-57. Because many of these sources are transient in nature, the 
continuous monitoring capability of BATSE allowed timing measurements to 
be performed whenever the systems were active, thus providing the necessary 
coverage to determine or improve the orbital parameters. 

Table 2 provides the current best-fit orbital parameters for 15 accreting 
pulsar systems. We have used BATSE data to derive orbital parameters for 9 
of these systems. We discuss a few of these briefly below. 


OAO 1657 - 415 . This system was initially discovered as an X-ray source by the 
Copernicus satellite (Polidan et al. 1978). The source was observed by several 
other X-ray satellites and showed both spin-up and spin-down behavior but no 
definite indication of orbital modulation of the pulse period. Observations with 
BATSE (Chakrabarty et al. 1993) detected the 10.4 day orbital period of this 
system and determined that OAO 1657-415 was an eclipsing high-mass system. 
The optical companion has not yet been identified, but the spectral type is likely 
B0-6Iab, inferred from the orbit and eclipse measurements. Identification of the 
companion and measurement of its orbital Doppler curve would constrain the 
mass of the neutron star in this system. 


GS 0834 - 430 . This system was initially discovered in an outburst during 
February, 1990 by the WATCH detectors on the Granat spacecraft (Lapshov 
et al. 1992). Pulsations at 12.3 s were discovered by Ginga in November, 1990 
and an accurate position determined (Aoki et al. 1992). A 114 day periodicity 
in the outbursts was later reported by Granat/WATCH, indicating a likely or- 
bital period (Lapshov et al. 1992). BATSE observations were carried out using 
both occultation analysis for unpulsed flux, and standard timing analysis for 
the pulsed component and has been reported in Wilson et al. (1994a). BATSE 
has observed seven outbursts of GS 0834-430 since GRO began observations, 
allowing measurement of the orbital parameters of the system (see Table 2). 

Accurate measurement of the orbital eccentricity is important, but is com- 
plicated by the uncertainty in the decoupling of the accretion and Doppler in- 
duced pulse-frequency changes. If the eccentricity is very close to zero, the 
companion is probably a 2-5 giant that has circularized the orbit due to 
dissipation in its envelope. On the other hand, if the eccentricity is small but 
finite (e.g. 0.1-0. 2) the companion is almost certainly not a giant, but might 
possibly be a B or perhaps a Be star. An optical/IR companion has not yet 
been identified. 



TABLE 2 

Orbital Parameters ok Accreting Pulsar Systems 
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Orbital epoch for longitude 309° ± 15°. Held fixed in fitting other elements. 


EXO 2030+375. This system was discovered by EXOSAT in May, 1985 (Par- 
mar et al. 1989). The orbital period was determined to be ~ 46 days, with some 
ambiguity in the orbital determination due to the finite extent of the EXOSAT 
observations. BATSE has observed EXO 2030+375 in a dozen consecutive or- 
bits, allowing a very precise orbit to be determined, particularly for orbital 
phases near periastron (see Stollberg et al. 1994 and Table 2). 

A 0535+26. Although this system has been observed in outburst many times 
and at many wavelengths (for a review, see Giovannelli and Graziati, 1992), the 
orbit of the system has never been accurately determined. Recent activity in 
1993 of the A 0535+26 system has allowed observations over three consecutive 
orbits with BATSE leading to the first definitive published orbital parameters 
(see Finger et al. 1994b and Table 2). 

GRO J1008-57. This 93.5 s pulsar was discovered by BATSE on 14 July 
1993, reached maximum intensity approximately 10 days later, and remained 
active for a total of about one month (Wilson et al. 1994c). The pulsar showed 
significant frequency evolution with spin-down and spin-up behavior roughly 
correlated with the hard X-ray flux increase and decrease respectively, similar 
to the Doppler-induced behavior observed in outbursts of EXO 2030+375. The 
spin evolution of the system, its transient nature, the fact that the luminosity 
is a significant fraction of the Eddington luminosity for a fiducial distance of 5 
kpc, and the low galactic latitude of 1° are all strongly suggestive of a Be-binary 
pulsar. From the Corbet-diagram (Figure 1) a 93.5 s Be-binary pulsar might 
be expected to have an orbital period of 100-200 d. However, no additional 
emission has been detected from this source through the end of 1993. 

Her X-l. This system is perhaps the best studied of all accreting binary pulsars. 
New results from BATSE (Wilson et al. 1994b) show a correlation between the 
hard X-ray flux at the peak of the main-on portion of the 35d cycle and the 
spin-frequency derivative, i.e. spin-down is correlated with low-flux as predicted 
by models such as those of Ghosh and Lamb ( 1979). Also, BATSE observations 
of the turn-on times of the 35d cycle indicate a possible correlation of early 
turn-on with decreased mass- transfer. 

4U 1538-52. Analysis of this system using BATSE data has only recently 
begun. Preliminary results given by Rubin et al. (1994) show a secular spin-up 
trend after almost 10 years of spin-down behavior. Further BATSE analysis will 
yield a significantly improved orbital solution for this system. 

Torque Studies 

Accreting pulsars provide a laboratory for the studies of the torque ex- 
hibited during magnetic accretion, allowing for a comparison to existing theory 
(see Ghosh & Lamb 1979). The continuous monitoring capability of BATSE has 
allowed long-term studies of accretion torques in several systems. These include 
Her X-l, Vela X-l, C’en X-3, GX 1+4, 4U 1626-67. and OAO 1657-415. The 
frequency histories for these sources for a period of over two years are shown 
in Figures 3a&b. Figure 3a shows spin-frequency histories for sources whose 
accretion is thought to be disk-fed. These are discussed in more detail below. 
For comparison, Figure 3b shows measurements of Vela X-l, which is a wind-fed 
accretor, and OAO 1657-415, which may be either disk-fed or wind-fed and sits 
in the Corbet diagram intermediate between the two types of systems. 
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Figure 3a. Spin-frequency history for probable disk-fed systems. 

(For Her X-l, mean frequencies for each main-on portion of the 35d cycle.) 



Figure 3b. Spin-frequency history for probable wind-fed systems. 









The very different qualitative behavior of the various systems is immedi- 
ately apparent from Figure 3. In particular, GX 1+4 and 4U i626-67 both show 
a monotonic spin-down with nearly constant average rate while other disk-fed 
and wind-fed systems exhibit frequent changes in sign of the frequency deriva- 
tive. 

The torque however does switch sign on decade time scales in GX 1+4 and 
4U 1626-67. From its discovery until the late 1980s, GX 1+4 was observed to 
be in a state of monotonic spin-up. Given its short spin-up time ( t su ~ 40 yr) 
compared to the system evolution time {t evo i > 10 6 yr), this was clearly not a 
permanent state, and GINGA and BATSE observations of a switch to spin-down 
in GX 1+4 (Makishima et al. 1988 Chakrabarty et al. 1994) showed that this 
was indeed the case. BATSE has found that 4U 1626-67 (t su ~ 4000 yr) has also 
changed from spin-up to spin-down ( Bildsten et al. 1994), and surprisingly, with 
the same average value of the torque but opposite sign. Since a disk reversal in 
this Roche-lobe overflow system is hard to imagine, we interpret the steady spin- 
down as a sign that the pulsar is spinning near its equilibrium period, where the 
magnetospheric radius equals the co- rotation radius (see Ghosh & Lamb 1979), 
requiring a magnetic field strength of order % 3 x 10 12 G for 4U 1626-67 (and 
10 14 G for GX 1+4). 

The simplest working hypotheses given the observations and the accretion 
torque theories of Ghosh ancl Lamb are then: (1) since t su <C tevoh X-ray 
pulsars quickly reach the equilibrium period, as determined from the accretion 
rate and magnetic field strength, (2) once there, the X-ray pulsar oscillates 
about the equilibrium period by having counteracting periods of spin-up and 
spin-down. This hypothesis was checked with pre-BATSE data by comparing 
the observed long-term (> vrs) torque to the fiducial value, Nf — M(GMr C0 ) 1 / 2 , 
where M is the maximum inferred mass accretion rate, M is the neutron star 
mass, and r co = (GM/u> 2 ) 1 ^ is the co-rotation radius and u> s is the angular 
spin frequency. For both GX 1+4 and 4L T 1626-67 (and also SMC X-l and IE 
2259+586) the observed long-term average torque, iV 0 , always satisfies N 0 > 
0.2 Nf, which, within the uncertainties, is consistent with accretion from matter 
near the co-rotation radius. These sources, which might be called ‘‘‘steady- 
staters”, show monotonic spin-up or spin-down over a few year time scale. 

As can be seen from Figure 3, other disk-fed sources such as Her X-l 
and Cen X-3 show quite different behavior. The long-term (> yrs) torque, 
N 0 , measured for these sources (as well as LMC X-4) was nearly a factor of 
100 smaller than Nf, implying different torquing behavior than in the “steady- 
staters”. These “wanderers” have frequency histories that are more consistent 
with a random- walk (Baykal & Ogelman 1993), which could arise if the torque 
had a value near Nf, but changed sign on a shorter timescale. The measured 
long term torque would thus be less than Nf. Recent short timescale (<10 — 20 d) 
torque measurements for Her X-l and Cen X-3 by BATSE always find a torque 
larger than N 0 , but never in excess of Nj (Wilson et al. 1994b, Finger et al. 
1994a). 

This classification scheme suggests that all disk-accreting pulsars show 
torques with magnitude <iVy on the shortest timescales and differentiate them- 
selves by the torque switching time. The wanderers switch within ~ 60 days, 
whereas the steady-staters switch once in 10-20 years. The primary issue for 
these systems is identifying the physics that sets this timescale. 
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